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A first order phase transition at high baryon density implies that a mixed phase can occupy 
a significant region of the interior of a neutron star. In this article we investigate the effect of a 
droplet phase on neutrino transport inside the core. Two specific scenarios of the phase transition 
are examined, one having a kaon condensate and the other having quark matter in the high density 
phase. The coherent scattering of neutrinos off the droplets greatly increases the neutrino opacity of 

. . . the mixed phase. We comment on how the existence of such a phase will affect a supernova neutrino 

0^ ' signal. 

a^ : 

0> ■ PACS numbers(s): 13.15.+g, 26.60.+C, 97.60.Jd 

O ^. Neutrino interactions play a central role in the early evolution of neutron stars formed during a type II supernova. 

D ' The enormous binding energy (~ 10^^ ergs) gained during implosion of the inner core of a massive star is released in 
C/^ , the form of neutrinos. Temporal and spectral characteristics of the neutrino emission from the newly born neutron 
f — ■ star, also called a protoneutron star, influences the explosion mechanism, the r-process nucleosynthesis, and more 
T-H , importantly, is an observable feature of a supernova explosion. The handful of neutrino events observed from SN87A 
are testimony to this fact. 

" ^ , In this article we study how a first order phase transition at high density will influence the neutrino opacity. We 
' consider two models of high density matter wherein a first order phase transition occurs for densities of relevance for 
^ neutron stars: first order kaon condensation and the quark-hadron transition. The mixed phase consists of a high 
baryon density, negatively charged, kaon condensed matter or quark matter coexisting with lower density, positively 
, charged, baryonic matter. We study the region of the mixed phase where the structure is expected to be droplets of 
kaonic matter or quark matter embedded in a lower baryon density nucleonic matter. Our primary motivation for 
this study is to determine the neutrino opacity of this material since this is directly related to an observable - the 
supernova neutrino luminosity curve. 

Our main finding is that the neutrino mean free path in a mixed phase containing either kaon condensed matter 
or quark matter is greatly reduced compared to uniform baryonic matter at the same density. This is a consequence 
O , of the coherent scattering of neutrinos from the matter in the droplets. The droplets carry a large weak charge and 
2 ' scattering of long-wavelength neutrinos increases as its square. The reduction of the mean free path compared to that 
in pure neutron matter could be as large as one to two orders of magnitude. Thus, despite our poor knowledge of 
neutrino opacities in dense nucleonic matter, we argue that the presence of a first order phase transition is likely to 
have a dramatic and discernible effect on the temporal characteristics of supernova neutrino light curve. 
; ^ Neutrino- Droplet Scattering: Prior to discussing in detail the specific models in which a droplet phase is energetically 

5^ , favored, we describe neutrino-droplet scattering. Neutrinos scatter off the droplets as they carry a net excess of weak- 
charge. For the typical droplet sizes which range from 5 — 15 fm we may infer that the scattering will be fairly coherent 
for momentum transfers ^ 40 MeV. 

First, we consider neutrino scattering from isolated droplets. The Lagrangian that describes the neutral current 
coupling of neutrinos to the droplet is given by 
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where Jj^ is the neutral current carried by the droplet and Gp = 1.166 x 10^*^ GeV^^ is the Fermi weak coupling 
constant. For non-relativistic droplets — p\y{x) S'^'"'^ has only a time like component and pwi^) is the excess weak 
charge density in the droplet. The total weak charge enclosed in a droplet of radius r^ is given by Nw — Jq'' d^x pw{x) 
and the form factor is F{q) = {1/Nw) Jq'' (Px pw{x) sinqx/qx. The differential cross section for neutrinos scattering 
from an isolated droplet is then given by 

_^ = iG|<(l + co».)i.'(„. (2) 

In the above equation E^, is the neutrino energy and 6 is the scattering angle. Since the droplets are massive we 
consider only elastic scattering in which case the magnitude of the momentum transfer is given hy q — \/2E^{l — cos9). 
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Neutrino transport in the mixed phase: We must also embed the droplets into a description of the medium to 
evaluate the neutrino transport parameters. The droplet radius ra and the inter-droplet spacing are determined by 
the interplay of surface and Coulomb energies and depends on the specific details of the high density model. We will 
calculate the droplets in the Wigner-Seitz cell assuming that the medium has a droplet density No- The cell radius 
Rw = (3/47riV£))^/'^. Except for one aspect, we will neglect the coherent scattering from more than one droplet. If 
the droplets form a lattice, Bragg scattering will dominate and our description would not be valid. But for a low 
density and a liquid phase the interference from multiple droplets affects only the scattering at long wavelengths. If 
the ambient temperature is not small compared to the melting temperature the droplet phase will be a liquid and for 
neutrino energy in the range (l/Rw) interference effects arising from scattering off different droplets is small. 
However, multiple droplet scattering cannot be neglected for 1/Rw The effects of other droplets is to cancel 
the scattering in the forward direction, because the interference is destructive except at exactly zero degrees, where it 
produces a change in the index of refraction of the medium. We account for this by subtracting from the weak charge 
density pw a uniform density which has the same total weak charge Nw as the matter in the Wigner-Seitz cell. The 
effect on the form factor is to modify it as follows 



F{q) ^ Fiq) = Fiq) - 3 



sin qRw — (qRw) cos qRw 



(3) 



In the medium, the appropriate differential cross section is defined per unit volume rather than per droplet. The 
neutrino-droplet differential cross section per unit volume follows from the preceding discussion and is given by 
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Note that even for small droplet density Njj the factor acts to enhance the droplet scattering. To quantify the 
importance of droplets as a source of opacity we compare with the standard scenario where matter is uniform and 
composed of neutrons. In this case the dominant source of opacity is due to scattering from thermal fluctuations and 
the cross section per unit volume is given by 
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where cv,ca are the vector and axial coupling constants of the neutron, n„ is neutron number density, Efn is the 
neutron Fermi energy and T is the matter temperature . 

The transport cross sections that are employed in studying the diffusive transport of neutrinos in the core are 
weighted by the angular factor (1 — cos 9). The transport mean free path \{Ey) for given neutrino energy E^ is given 

by 
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It is thus the task of the models to provide the weak charge and form factors of the droplets so that we may calculate 
the contribution of neutrino-droplet scattering to the opacity of the mixed phase. 

First-order Kaon condensation Kaon condensation is one of the several possible forms of exotica that could exist 
at high baryon density. Since the pioneering work of Kaplan and Nelson |^ several authors (see and references 
therein) have studied in detail the role of such a condensate in neutron star structure and evolution. Kaon condensation 
significantly softens the high density equation of state (EOS) thereby lowering the maximum mass of neutron stars. 
Glendenning and Schaffner revisited the problem and studied kaon condensation in a relativistic mean field model, 
wherein the kaon is minimally coupled to a, lo and p mesons. Within the purview of this model they were able to 
show that kaon condensation occurs as a first order phase transition. 

A relativistic field theoretical approach is used to describe the baryons. We follow the model employed by Glen- 
denning and Schaffner |^ wherein the Lagrangian density is given by 

-Z^ = ^ 5(^7^9^ - 5c^B7°^o - gpB7°bo • t 
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Here, B are the Dirac field operators for baryons and t is the isospin matrix. In the mean field approximation only the 
zeroth components of the sigma, omega and rho meson fields are retained and are denoted by a, uq and bo respectively. 
U{(t) represents the scalar self-interactions and is of the form U{(t) — {b/3)Mn{gaN<^)^ + [c/ 4)[g„N(jY . Kaons are 
coupled minimally to the meson fields appearing in the baryonic part of the Lagrangian. The Lagrange density for 
the kaon field is given by 

Lk = VIK*V^K - (rriK - gaK<TfK*K , 

where the covariant derivative 



= [d^ + iigujK^^fj. + gpRbfj.)) ■ 

Droplets are constructed in the Thomas-Fermi or local density approximation. The equation of motion for the kaons 
is given by 

[D^V^ + {ml)^]K~ =0, (7) 

where m|f = rufc — gaK<^ is the kaon effective mass in the medium. Eq. is solved by writing explicitly the 
spatial and temporal dependence of the kaon field as K~ = /^0(r) exp(i^xt)/-\/2 where = 93 MeV is the pion 
decay constant. The spatial part has only a radial dependence, and the time dependence is fixed by the chemical 
potential 0. In principle we should solve a Klein-Gordon equation for the other meson fields as well. However, the 
other mesons are significantly more massive and their fields do not vary over distances of the order of their Compton 
wavelengths. In the vicinity of the critical density the effective kaon mass is approximately given by /Zg ~ 200 

MeV and is small compared to a, lo ot p meson masses. In this approximation the non-strange meson field equations 



are 



mlcT = + g^Bipn + Pp) + Q'ykf^O^ml , (8) 



-gUPn + Pp)~glKfn(}\m + X), (9) 



'T^pbo = T;9piPp " Pn) - gpKflO'^iPK + X) . (10) 



where pn, pp, p^ and p^ are the baryon number and scalar densities of the neutrons and protons respectively and the 
vector potential for the kaons is given hy X = g^jk^o + gpkbo- 

The classical kaon field equation is obtained by substituting the specific form given above in Eq. (Q) and is given by 

- V^0{r) + [{pK + Xf - (TOfc - gaaofW) = . (11) 

The above differential equation is solved with appropriate boundary conditions to obtain the droplet configuration. 
The baryon and electron chemical potentials are close to those in the mixed phase, but must be fine tuned to give 
a solution to Eq.(pT|) that satisfies the boundary conditions at both r = and r — Rw. The energy is computed in 
the Wigner-Seitz approximation and the droplet radius that minimizes the energy density is determined. The energy 
density contains bulk, surface and Coulomb contributions. The baryonic contribution is 

eB [r) = efe,„ + ^mlaj + ^m^ ^ 1 ^2^2 + [/(^^) ^ (12) 

and the kaon contribution is given by 

ek{r)=f^e'ir){t,K+Xf. (13) 

Note that the surface contribution is contained in the above expression and that terms proportional to the gradients 
have been eliminated by using the equation of motion for the kaons. 

The shape of the cell depends on the type of lattice that is favored; for BCC or FCC the Wigner-Seitz cell is a 
regular polyhedron and can be approximated by a sphere. The Wigner-Seitz radius, denoted by Rw, is determined 
by size of the cell that encloses zero total electric charge. Since the electric field vanished at r = Rw, the Coulomb 
energy is given by 

P a /-^'^ , Qlir) ,^ 
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where a = 1/137 and Qeif) is the enclosed electric charge a distance r from the center. 

We find that the Debye screening lengths (1/A|))i = 4:Tra{dni/dfii)nj,j^i of electrons, protons and kaons are of the 
same order as that the typical droplet radius. Thus the charged particle distributions may change and screen the 
droplet charge We allow for this by adding the Coulomb term to the single particle energies and employing a 
iterative procedure to relax the charged particle distributions to their equilibrium value. For droplets radii rd ^ \d 
the charge particle distributions in the droplet changes at the 10-30 percent level depending on their individual Debye 
screening lengths. However, for rd/Xo ^ 1 the charge particle distributions change drastically and all the charge 
resides in the surface region. For these conditions the droplet phase in not energetically favored . Our investigation 
indicates that for baryon density in the region 0.52 — 0.62 fm~'^ the droplets survive, but with moderate modifications 
to the charged particle profiles. Fig. 1 shows the particle density profiles at ub — 0.54 fm~^. 
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FIG. 1. Droplet profile at baryon density ub = 0.54 fm~^ . The solid curve shows the kaon number density, the short-dashed 
and long-dashed curves show the proton and neutron number densities respectively. The kaon, proton, neutron and electron 
number densities are shown as solid, dashed, long-dashed and dot-dashed curves respectively. 



The excess vector weak charge density of the droplet is 

Pw{r) = {-1 + 2 sin^ 9w)nK{r) 

+ (l + 4sin2 0w)ine{r) ~ n°) 

+ (1 - 4sin2 ew)inp{r) - n°) - K(r) - <) , (15) 

where uk-, Ue, Up and n„ are the kaon, electron, proton and neutron number densities respectively and n°, n° and n° 
are corresponding densities in the low density phase. Kaons dominate the weak charge density since (i) the electron 
number is small, (ii) the protons carry negligible weak charge, and (iii) the neutron number density inside and outside 
the droplet are nearly equal. The neutral current neutrino-kaon coupling is not experimentally measured and our 
estimates here are based on the constituent quark model (note that the main source of uncertainty is the singlet 
contribution in the large Nc limit that the singlet current contributes at the 30 percent level to the weak charge 
of the kaon). The total weak charge of the droplet shown in Fig. 1 is Nw — 700. The corresponding form factor F{q) 
is shown in the top panel of Fig. 2. 

The transport differential scattering (Eq.(||)) weighted with the factor (1 — cos6') for E^, — 50 MeV and for a baryon 
density of ns — 0.54 fm^'^ is shown in the bottom panel of Fig. 2. For comparison we also show results for uniform 
neutron matter (computed by using Eq.(^ and for T = 10 MeV) where thermal density fluctuations dominate the 
opacity. 
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FIG. 2. Top panel: The dashed curve shows the form factor for a single drop and the solid curve shows the results for an 
embedded droplet (Eq.(^). Lower panel: Transport differential cross section for scattering off droplets (solid curve) and for 
scattering off thermal fluctuations in pure neutron matter (dashed curve). Note that dashed curve has been multiplied by a 
factor of 10. 

The neutrino mean free path computed using Eq.(^) is shown (solid curve) in Fig. 3 for ub — 0.54 and for a 
temperature of 10 MeV. The mean free path in uniform neutron matter is also shown (dashed curve) for comparison. 
As is evident from the figure, the mean free path in the mixed phase is significantly reduced compared to pure 
neutron matter for typical neutrino energy Ei, ^ ttT. The reduction is severe at moderate energies, and is easily 
understood by noting that the form factors have significant support in the region where g ~ 40 MeV. At lower q the 
inter-droplet correlations act to screen the weak charge of the droplet and at higher energy is attenuated by the form 
factor. Although the transport calculations would require the energy dependent mean free path it is useful to define 
an average mean free path for energy transport < A^; >= J dE,y fu{E^) Ef;X{E^)/{J dE^ fi,{Ei,) E^) where fu{E^) is 
the neutrino distribution function At T = lOMel^we find that < A^; > is reduced by a factor ~ 20 and at T = 30 
MeV the reduction is reduced by a factor ~ 9. 
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FIG. 3. Neutrino mean free path in the droplet phase. Results for neutrino droplet scattering (solid curve) are compared 
with the uniform neutron matter case (dashed curve). 
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First order quark-hadron transition Several authors have studied the possibihty of a first order quark-hadron phase 
transition in a bag model picture for densities of relevance to neutron star interiors [pl]-p^ . These studies found that 
a mixed phase was energetically favored. Detailed investigations of the structured phase by Heiselberg et ai, wherein 
surface and Coulomb contributions to the energy density were properly accounted for, showed that the presence of 
a droplet phase was sensitive to the value of the surface tension a between the quark and hadron phases Q . For 
a < 100 MeV the droplet sizes favored were smaller than the typical Debye screening length Xd ^ 7 fm and the 
droplet phase was energetically favored. 

1/3 

The droplet radius is given by = [l5cr/(87r(pQ — pn)'^)] , where pg and are the electric charge densities 
in the quark and nucleon phases respectively. The charge densities depend in general on details of the model and a 
typical scale for the difference pg — pN ~ 0.4 e fm~^ . For a in the range 10 — 100 MeV this corresponds to droplet radii 
of = 3 — 7 fm. For our estimate of coherent scattering we will assume a droplet radius = 5 fm at a baryon density 
ub = 0.7 fm~'^. In a minimal model the nuclconic phase is described by a relativistic mean field theory and the quark 
phase as a free Fermi gas with a bag constant B — 200 MeV/fm^ and a finite strange quark mass rUg = 150 MeV. 
In this model, and for ub = 0.7 fm~^, the baryon density in the nucleon phase is njv = 0.68 fm~^ and uq — 0.94 fm~'^ 
in the quark phase. The volume fraction occupied by the quark phase is given by / = (jib — nN)/{nQ — njv) ~ 0.1 
and the Wigner-Seitz cell radius spacing Rw = r^l j^l'^ ~ 11 fm. 

For small droplets the quark densities are uniform to good approximation and are determined by conditions of local 
chemical equilibrium. The electron density is the same inside and outside the droplet. Thus the excess weak charge 
density in the quark droplet is given by 

Pw = [("« - + n^)) - (n„ - n^)\ 

^ ^["u ^ ("d + "-s) + sin^ 6*1^ , (16) 

where nu,nd,ns are the up, down and strange quark number densities inside the droplet, and Un^Up the neutron and 
proton densities outside. The individual particle densities are shown in the Fig. ^. The net weak charge and form 
factor for a droplet of radius = 5 fm and cell radius Rw = 11 fm are computed as described earlier. We find that 
the total weak charge Nw ~ 860 and the droplet density No = 18 x 10^^ fm^'^. The transport differential cross 
section (left panel) and the neutrino mean free path (right panel) are shown in Fig. 5. We see that the results are 
qualitatively similar to those in the condensate scenario. The energy averaged neutrino mean free path < A^; > 
defined earlier at T = 10 MeV is reduced by a factor of ~ 10 compared to pure neutron matter and at T = 30 MeV 
is reduced by a factor ^ 17. 
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FIG. 4. Results for neutrino-quark droplet scattering (solid curve) are compared with the uniform neutron matter case 
(dashed curve). Left panel shows the results for the differential scattering of neutrinos with energy Ev = 50 MeV and the right 
panel shows the neutrino mean free path at ns ~ 0.7 fm"'^ and T = 10 MeV. The left panel also shows the individual particle 
densities in the mixed phase. 
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Discussion: There are several caveats to the study presented in this work which we would like to address in this 
paragraph. First, formula for the scattering in Eq.(||) treats the effects of droplet-droplet correlations rather crudely, 
and short range order between droplets could affect the low momentum transfer scattering. However, these momenta 
have small weight at temperatures of interest, and corrections to the transport would be small. Our treatment of the 
form factor at small q may be viewed as being conservative. Second, in the mixed phase we studied only droplets 
and ignored the possibility of other spatial structures which might exist for some range of density. Qualitatively, the 
heterogeneity of the mixed phase is the key feature that results in an enhanced neutrino opacity. Thus, our results 
provide a rough estimate of what one may expect even when the structures are more complex. Both these issues 
clearly require more attention and we hope to pursue this in detail in future work. Thirdly, the droplets are likely to 
have low lying collective modes which the neutrinos could excite. This will also make an important contribution to 
the opacity of the material and requires further work. Lastly, we expect that droplet formation will be fairly rapid 
due to the large temperatures (T ~ 50 MeV) reached during the implosion of the inner core. 

Our main focus was to study the consequences of a first order phase transition on the transport of neutrinos in the 
inner core of neutron stars. We have shown that the droplet phase is indeed very opaque to neutrinos. Studies of 
the early evolution of neutron stars, where diffusive transport of neutrinos directly affects the temporal evolution of 
neutrino emission, which incorporate the effects arising due to a first order phase transition remain largely unexplored. 
In view of our findings in this article we may expect significant changes to the neutrino emission if a droplet phase were 
to exist in the inner core of a newly born neutron star. It seems likely that detailed studies of the diffusive transport of 
neutrinos in models with a droplet phase will provide the link between the observable features of a supernova neutrino 
signal and the phases of matter at high density. 

The authors would like to thank Paulo Bedaque, Chuck Horowitz, David Kaplan, Thomas Papenbrock, Martin 
Savage and Ray Sawyer for useful discussions. We would also like to thank Jason Cooke for providing us numerical 
tables of his work on quark-hadron phase transitions. This work is supported by the grant DOE-FG-06-90ER40561. 
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